Discovery of a new photosynthetic euglenoid in Poland: Euglena mazurica sp. nov. (Euglenales, Euglenaceae).

Herein we describe a new photosynthetic euglenoid species found in Poland - Euglena mazurica. A large population exists in a small, eutrophic body of water located in a pasture near Mikolajki town inside the Masurian Landscape Park (covering a part of the Masurian Lake District in Poland). The unique cell shape (corkscrew-like) discerns it well from other previously described euglenoid species with metabolic cells. The new species possesses two plate-like chloroplasts each with a pyrenoid accompanied by two paramylon caps placed on either side of it (diplopyrenoids). On the phylogenetic tree, the new species is situated within the Euglena clade. Though it is a sister branch to three clades - one representing the similar Euglena agilis, characterized by its fusiform cells and two chloroplasts with diplopyrenoids, the two species are clearly morphologically and molecularly distinct.

Single-cell genomics revealed Candidatus Grellia alia sp. nov. as an endosymbiont of Eutreptiella sp. (Euglenophyceae).

Though endosymbioses between protists and prokaryotes are widespread, certain host lineages have received disproportionate attention what may indicate either a predisposition to such interactions or limited studies on certain protist groups due to lack of cultures. The euglenids represent one such group in spite of microscopic observations showing intracellular bacteria in some strains. Here, we perform a comprehensive molecular analysis of a previously identified endosymbiont in the Eutreptiella sp. CCMP3347 using a single cell approach and bulk culture sequencing. The genome reconstruction of this endosymbiont allowed the description of a new endosymbiont Candidatus Grellia alia sp. nov. from the family Midichloriaceae. Comparative genomics revealed a remarkably complete conjugative type IV secretion system present in three copies on the plasmid sequences of the studied endosymbiont, a feature missing in the closely related Grellia incantans. This study addresses the challenge of limited host cultures with endosymbionts by showing that the genomes of endosymbionts reconstructed from single host cells have the completeness and contiguity that matches or exceeds those coming from bulk cultures. This paves the way for further studies of endosymbionts in euglenids and other protist groups. The research also provides the opportunity to study the diversity of endosymbionts in natural populations.

Gleaning Euglenozoa-specific DNA polymerases in public single-cell transcriptome data.

Multiple genes encoding family A DNA polymerases (famA DNAPs), which are evolutionary relatives of DNA polymerase I (PolI) in bacteria and phages, have been found in eukaryotic genomes, and many of these proteins are used mainly in organelles. Among members of the phylum Euglenozoa, distinct types of famA DNAP, PolIA, PolIBCD+, POP, and eugPolA, have been found. It is intriguing how the suite of famA DNAPs had been established during the evolution of Euglenozoa, but the DNAP data have not been sampled from the taxa that sufficiently represent the diversity of this phylum. In particular, little sequence data were available for basal branching species in Euglenozoa until recently. Thanks to the single-cell transcriptome data from symbiontids and phagotrophic euglenids, we have an opportunity to cover the "hole" in the repertory of famA DNAPs in the deep branches in Euglenozoa. The current study identified 16 new famA DNAP sequences in the transcriptome data from 33 phagotrophic euglenids and two symbiontids, respectively. Based on the new famA DNAP sequences, the updated diversity and evolution of famA DNAPs in Euglenozoa are discussed.

Molecular phylogenetics of sessile Dolium sedentarium, a petalomonad euglenid.

The euglenids are a species-rich group of flagellates with varying modes of nutrition that can be found in diverse habitats. Phagotrophic members of this group gave rise to phototrophs and hold the key to understanding the evolution of euglenids as a whole, including the evolution of complex morphological characters like the euglenid pellicle. Yet to understand the evolution of these characters, a comprehensive sampling of molecular data is needed to correlate morphological and molecular data, and to estimate a basic phylogenetic backbone of the group. While the availability of SSU rDNA and, more recently, multigene data from phagotrophic euglenids has improved, several "orphan" taxa remain without any molecular data whatsoever. Dolium sedentarium is one such taxon: It is a rarely-observed phagotrophic euglenid that inhabits tropical benthic environments and is one of few known sessile euglenids. Based on morphological characters, it has been thought of as part of the earliest branch of euglenids, the Petalomonadida. We report the first molecular sequencing data for Dolium using single-cell transcriptomics, adding another small piece in the puzzle of euglenid evolution. Both SSU rDNA and multigene phylogenies confirm it as a solitary branch within Petalomonadida.

Taxonomy of a New Parasitic Euglenid, Euglenaformis parasitica sp. nov. (Euglenales, Euglenaceae) in Ostracods and Rhabdocoels.

Parasitic euglenids have rarely been studied. We found parasitic euglenids in two species of ostracods (Cyprinotus cassidula, Dolerocypris sinensis) and two species of rhabdocoels (Mesostoma lingua, Microdalyellia armigera) in a rice field. These parasites grew and proliferated inside the host body. These parasites had pellicle strips, one emergent flagellum, and a red stigma, but no chloroplasts, and showed euglenoid movement. Inside the living host, they did not have emergent flagella and moved only by euglenoid movement, but when the host died or the parasites were isolated from the host, they extended their flagella and switched to swimming movement. We conclude that the parasites found in the four hosts that we examined are of the same species, considering the morphological characteristics and identities in the nSSU and nLSU rDNA sequences of those parasites. Molecular phylogenetic analysis showed that the parasite formed a clade with the free-living photoautotrophic species of Euglenaformis, with moderate statistical support. Therefore, the parasite is a secondary osmotroph derived from a photoautotrophic ancestor. Based on the results of morphological observation and molecular phylogenetic analysis, we propose a new species of parasitic euglenid, Euglenaformis parasitica sp. nov.

Hydrophobicity-Driven Increases in Editing in Mitochondrial mRNAs during the Evolution of Kinetoplastids.

Kinetoplastids are a diverse group of flagellates which exhibit editing by insertion/deletion of Us in the mitochondrial mRNAs. Some mRNAs require editing to build most of their coding sequences, a process known as pan-editing. Evidence suggests that pan-editing is an ancestral feature in kinetoplastids. Here, we investigate how the transition from nonedited to pan-edited states occurred. The mitochondrial mRNAs and protein sequences from nine kinetoplastids and related groups (diplonemids, euglenids, and jakobids) were analyzed. RNA editing increased protein hydrophobicity to extreme values by introducing Us in the second codon position, despite the absence of editing preferences related to codon position. In addition, hydrophobicity was maintained by purifying selection in species that lost editing by retroposition of the fully edited mRNA. Only a few hydrophobic to hydrophilic amino acid changes were inferred for such species. In the protein secondary structure, these changes occurred spatially close to other hydrophilic residues. The analysis of coevolving sites showed that multiple changes are required together for hydrophobicity to be lost, which suggest the proteins are locked into extended hydrophobicity. Finally, an analysis of the NAD7 protein-protein interactions showed they can also influence hydrophobicity increase in the protein and where editing can occur in the mRNA. In conclusion, our results suggest that protein hydrophobicity has influenced editing site selection and how editing expanded in mRNAs. In effect, the hydrophobicity increase was entrenched by a neutral ratchet moved by a mutational pressure to introduce Us, thus helping to explain both RNA editing increase and, possibly, persistence.

Phylogenomics of novel ploeotid taxa contribute to the backbone of the euglenid tree.

Euglenids are a diverse group of flagellates that inhabit most environments and exhibit many different nutritional modes. The most prominent euglenids are phototrophs, but phagotrophs constitute the majority of phylogenetic diversity of euglenids. They are pivotal to our understanding of euglenid evolution, yet we are only starting to understand relationships amongst phagotrophs, with the backbone of the tree being most elusive. Ploeotids make up most of this backbone diversity-yet despite their morphological similarities, SSU rDNA analyses and multigene analyses show that they are non-monophyletic. As more ploeotid diversity is sampled, known taxa have coalesced into some subgroups (e.g. Alistosa), but the relationships amongst these are not always supported and some taxa remain unsampled for multigene phylogenetics. Here, we used light microscopy and single-cell transcriptomics to characterize five ploeotid euglenids and place them into a multigene phylogenetic framework. Our analyses place Decastava in Alistosa; while Hemiolia branches with Liburna, establishing the novel clade Karavia. We describe Hemiolia limna, a freshwater-dwelling species in an otherwise marine clade. Intriguingly, two undescribed ploeotids are found to occupy pivotal positions in the tree: Chelandium granulatum nov. gen. nov. sp. branches as sister to Olkasia, and Gaulosia striata nov. gen. nov. sp. remains an orphan taxon.

18S rDNA sequence-structure phylogeny of the Euglenophyceae (Euglenozoa, Euglenida).

The phylogeny of Euglenophyceae (Euglenozoa, Euglenida) has been discussed for decades with new genera being described in the last few years. In this study, we reconstruct a phylogeny using 18S rDNA sequence and structural data simultaneously. Using homology modeling, individual secondary structures were predicted. Sequence-structure data are encoded and automatically aligned. Here, we present a sequence-structure neighbor-joining tree of more than 300 taxa classified as Euglenophyceae. Profile neighbor-joining was used to resolve the basal branching pattern. Neighbor-joining, maximum parsimony, and maximum likelihood analyses were performed using sequence-structure information for manually chosen subsets. All analyses supported the monophyly of Eutreptiella, Discoplastis, Lepocinclis, Strombomonas, Cryptoglena, Monomorphina, Euglenaria, and Colacium. Well-supported topologies were generally consistent with previous studies using a combined dataset of genetic markers. Our study supports the simultaneous use of sequence and structural data to reconstruct more accurate and robust trees. The average bootstrap value is significantly higher than the average bootstrap value obtained from sequence-only analyses, which is promising for resolving relationships between more closely related taxa.

Typical structure of rRNA coding genes in diplonemids points to two independent origins of the bizarre rDNA structures of euglenozoans.

BACKGROUND: Members of Euglenozoa (Discoba) are known for unorthodox rDNA organization. In Euglenida rDNA is located on extrachromosomal circular DNA. In Kinetoplastea and Euglenida the core of the large ribosomal subunit, typically formed by the 28S rRNA, consists of several smaller rRNAs. They are the result of the presence of additional internal transcribed spacers (ITSs) in the rDNA. Diplonemea is the third of the main groups of Euglenozoa and its members are known to be among the most abundant and diverse protists in the oceans. Despite that, the rRNA of only one diplonemid species, Diplonema papillatum, has been examined so far and found to exhibit continuous 28S rRNA. Currently, the rDNA organization has not been researched for any diplonemid. Herein we investigate the structure of rRNA genes in classical (Diplonemidae) and deep-sea diplonemids (Eupelagonemidae), representing the majority of known diplonemid diversity. The results fill the gap in knowledge about diplonemid rDNA and allow better understanding of the evolution of the fragmented structure of the rDNA in Euglenozoa. RESULTS: We used available genomic (culture and single-cell) sequencing data to assemble complete or almost complete rRNA operons for three classical and six deep-sea diplonemids. The rDNA sequences acquired for several euglenids and kinetoplastids were used to provide the background for the analysis. In all nine diplonemids, 28S rRNA seems to be contiguous, with no additional ITSs detected. Similarly, no additional ITSs were detected in basal prokinetoplastids. However, we identified five additional ITSs in the 28S rRNA of all analysed metakinetoplastids, and up to twelve in euglenids. Only three of these share positions, and they cannot be traced back to their common ancestor. CONCLUSIONS: Presented results indicate that independent origin of additional ITSs in euglenids and kinetoplastids seems to be the most likely. The reason for such unmatched fragmentation remains unknown, but for some reason euglenozoan ribosomes appear to be prone to 28S rRNA fragmentation.

Maturyoshka: A maturase inside a maturase, and other peculiarities of the novel chloroplast genomes of marine euglenophytes.

Organellar genomes often carry group II introns, which occasionally encode proteins called maturases that are important for splicing. The number of introns varies substantially among various organellar genomes, and bursts of introns have been observed in multiple eukaryotic lineages, including euglenophytes, with more than 100 introns in their plastid genomes. To examine the evolutionary diversity and history of maturases, an essential gene family among euglenophytes, we searched for their homologs in newly sequenced and published plastid genomes representing all major euglenophyte lineages. We found that maturase content in plastid genomes has a patchy distribution, with a maximum of eight of them present in Eutreptiella eupharyngea. The most basal lineages of euglenophytes, Eutreptiales, share the highest number of maturases, but the lowest number of introns. We also identified a peculiar convoluted structure of a gene located in an intron, in a gene within an intron, within yet another gene, present in some Eutreptiales. Further investigation of functional domains of identified maturases show that most of them lost at least one of the functional domains, which implies that the patchy maturase distribution is due to frequent inactivation and eventual loss over time. Finally, we identified the diversified evolutionary origin of analysed maturases, which were acquired along with the green algal plastid or horizontally transferred. These findings indicate that euglenophytes' plastid maturases have experienced a surprisingly dynamic history due to gains from diversified donors, their retention, and loss.

Toward the robust resolution of taxonomic ambiguity within Lepocinclis (Euglenida) based on DNA sequencing and morphology.

DNA sequences were analyzed for three groups of species from the Lepocinclis genus (L. acus-like, L. oxyuris-like, and L. tripteris-like) along with cellular morphology. Phylogenetic analyses were based on nuclear SSU rDNA, LSU rDNA, and plastid-encoded LSU rDNA. DNA sequences were obtained from species available in culture collections (L. acus SAG 1224-1a and UTEX 1316) and those isolated directly from the environment in Poland (48 isolates), resulting in 79 new sequences. The obtained phylogenetic tree of Lepocinclis included 27 taxa, five of which are presented for the first time (L. convoluta, L. gracillimoides, L. longissima, L. pseudospiroides, and L. torta) and nine taxonomically verified and described. Based on morphology, literature data, and phylogenetic analyses, the following species were distinguished: in the L. acus-like group, L. longissima and L. acus; in the L. tripteris-like group, L. pseudospiroides, L. torta, and L. tripteris; in the L. oxyuris-like group, L. gracillimoides, L. oxyuris var. oxyuris, and L. oxyuris var. helicoidea. For all verified species, diagnostic descriptions were emended, nomenclatural adjustments were made, and epitypes were designated.

Measuring Photophysiology of Attached Stage of Colacium sp. by a Cuvette-Type Fast Repetition Rate Fluorometer.

Fast repetition rate fluorometer (FRRf) is a beneficial method for measuring photosystem II (PSII) photophysiology and primary productivity. Although FRRf can measure PSII absorption cross-section (σPSII), maximum photochemical efficiency (Fv/Fm), effective photochemical efficiency (Fq'/Fm'), and non-photochemical quenching (NPQNSV) for various eukaryotic algae and cyanobacteria, almost all FRRf studies to date have focused on phytoplankton. Here, the protocol describes how to measure PSII photophysiology of an epizoic alga Colacium sp. Ehrenberg 1834 (Euglenophyta), in its attached stage (attached to zooplankton), using cuvette-type FRRf. First, we estimated the effects of substrate zooplankton (Scapholeberis mucronata O.F. Müller 1776, Cladocera, Daphniidae) on baseline fluorescence and σPSII, Fv/Fm, Fq'/Fm', and NPQNSV of planktonic Colacium sp. To validate this methodology, we recorded photophysiology measurements of attached Colacium sp. on S. mucronata and compared these results with its planktonic stage. Representative results showed how the protocol could determine the effects of calcium (Ca) and manganese (Mn) on Colacium sp. photophysiology and identify the various effects of Mn enrichment between attached and planktonic stages. Finally, we discuss the adaptability of this protocol to other periphytic algae.

Heterotrophic euglenid Rhabdomonas costata resembles its phototrophic relatives in many aspects of molecular and cell biology.

Euglenids represent a group of protists with diverse modes of feeding. To date, only a partial genomic sequence of Euglena gracilis and transcriptomes of several phototrophic and secondarily osmotrophic species are available, while primarily heterotrophic euglenids are seriously undersampled. In this work, we begin to fill this gap by presenting genomic and transcriptomic drafts of a primary osmotroph, Rhabdomonas costata. The current genomic assembly length of 100 Mbp is 14× smaller than that of E. gracilis. Despite being too fragmented for comprehensive gene prediction it provided fragments of the mitochondrial genome and comparison of the transcriptomic and genomic data revealed features of its introns, including several candidates for nonconventional types. A set of 39,456 putative R. costata proteins was predicted from the transcriptome. Annotation of the mitochondrial core metabolism provides the first data on the facultatively anaerobic mitochondrion of R. costata, which in most respects resembles the mitochondrion of E. gracilis with a certain level of streamlining. R. costata can synthetise thiamine by enzymes of heterogenous provenances and haem by a mitochondrial-cytoplasmic C4 pathway with enzymes orthologous to those found in E. gracilis. The low percentage of green algae-affiliated genes supports the ancestrally osmotrophic status of this species.

Methods of Lipid Analyses for Microalgae: Charophytes, Eustigmatophytes, and Euglenophytes.

Algae are ecologically important organisms and are widely used for basic research, with a focus on for example photosynthesis, evolution, and lipid metabolism. Many biosynthetic pathways of algal lipids have been deciphered using available genomic information. Here we describe methods for lipid analyses from three representative algae, including Archaeplastida, the SAR lineage (Stramenopiles, Alveolata, Rhizaria), and Excavata. Archaeplastida acquired their plastids by primary endosymbiosis, and the others by secondary endosymbiosis with a Rhodophyceae-type plastid in SAR and a Chlorophyceae-type plastid in Excavata (Euglenozoa). Analytical methods for these algae are described for membrane lipids and neutral lipids including triacylglycerol and wax esters.

Single-cell genomics unveils a canonical origin of the diverse mitochondrial genomes of euglenozoans.

BACKGROUND: The supergroup Euglenozoa unites heterotrophic flagellates from three major clades, kinetoplastids, diplonemids, and euglenids, each of which exhibits extremely divergent mitochondrial characteristics. Mitochondrial genomes (mtDNAs) of euglenids comprise multiple linear chromosomes carrying single genes, whereas mitochondrial chromosomes are circular non-catenated in diplonemids, but circular and catenated in kinetoplastids. In diplonemids and kinetoplastids, mitochondrial mRNAs require extensive and diverse editing and/or trans-splicing to produce mature transcripts. All known euglenozoan mtDNAs exhibit extremely short mitochondrial small (rns) and large (rnl) subunit rRNA genes, and absence of tRNA genes. How these features evolved from an ancestral bacteria-like circular mitochondrial genome remains unanswered. RESULTS: We sequenced and assembled 20 euglenozoan single-cell amplified genomes (SAGs). In our phylogenetic and phylogenomic analyses, three SAGs were placed within kinetoplastids, 14 within diplonemids, one (EU2) within euglenids, and two SAGs with nearly identical small subunit rRNA gene (18S) sequences (EU17/18) branched as either a basal lineage of euglenids, or as a sister to all euglenozoans. Near-complete mitochondrial genomes were identified in EU2 and EU17/18. Surprisingly, both EU2 and EU17/18 mitochondrial contigs contained multiple genes and one tRNA gene. Furthermore, EU17/18 mtDNA possessed several features unique among euglenozoans including full-length rns and rnl genes, six mitoribosomal genes, and nad11, all likely on a single chromosome. CONCLUSIONS: Our data strongly suggest that EU17/18 is an early-branching euglenozoan with numerous ancestral mitochondrial features. Collectively these data contribute to untangling the early evolution of euglenozoan mitochondria.

Evaluation of V2 18S rDNA barcode marker and assessment of sample collection and DNA extraction methods for metabarcoding of autotrophic euglenids.

Even though the interest in metabarcoding in environmental studies is growing, euglenids are still underrepresented in both sea and freshwater bodies researches. The reason for this situation could be the unsuitability of universal eukaryotic DNA barcodes and primers as well as the lack of a verified protocol, suitable to assess euglenid diversity. In this study, using specific primers for the V2 hypervariable region of 18S rDNA for metabarcoding resulted in obtaining a high fraction (85%) of euglenid reads and species-level identification of almost 90% of them. Fifty species were detected by the metabarcoding method, including almost all species observed using a light microscope. We investigated three biomass harvesting methods (filtering, centrifugation and scraping the side of a collection vessel) and determined that centrifugation and filtration outperformed scrapes, but the choice between them is not crucial for the reliability of the analysis. In addition, eight DNA extraction methods were evaluated. We compared five commercially available DNA isolation kits, two CTAB-based protocols and a chelating resin. For this purpose, the efficiency of extraction, quality of obtained DNA, preparation time and generated costs were taken into consideration. After examination of the aforementioned criteria, we chose the GeneMATRIX Soil DNA Purification Kit as the most suitable for DNA isolation.

Multigene phylogenetics of euglenids based on single-cell transcriptomics of diverse phagotrophs.

Euglenids are a well-known group of single-celled eukaryotes, with phototrophic, osmotrophic and phagotrophic members. Phagotrophs represent most of the phylogenetic diversity of euglenids, and gave rise to the phototrophs and osmotrophs, but their evolutionary relationships are poorly understood. Symbiontids, in contrast, are anaerobes that are alternatively inferred to be derived euglenids, or a separate euglenozoan group. Most phylogenetic studies of euglenids have examined the SSU rDNA only, which is often highly divergent. Also, many phagotrophic euglenids (and symbiontids) are uncultured, restricting collection of other molecular data. We generated transcriptome data for 28 taxa, mostly using a single-cell approach, and conducted the first multigene phylogenetic analyses of euglenids to include phagotrophs and symbiontids. Euglenids are recovered as monophyletic, with symbiontids forming an independent branch within Euglenozoa. Spirocuta, the clade of flexible euglenids that contains both the phototrophs (Euglenophyceae) and osmotrophs (Aphagea), is robustly resolved, with the ploeotid Olkasia as its sister group, forming the new taxon Olkaspira. Ploeotids are paraphyletic, although Ploeotiidae (represented by Ploeotia spp.), Lentomonas, and Keelungia form a robust clade (new taxon Alistosa). Petalomonadida branches robustly as sister to other euglenids in outgroup-rooted analyses. Within Spirocuta, Euglenophyceae is a robust clade that includes Rapaza, and Anisonemia is a well-supported monophyletic group containing Anisonemidae (Anisonema and Dinema spp.), 'Heteronema II' (represented by H. vittatum), and a clade of Neometanema plus Aphagea. Among 'peranemid' phagotrophs, Chasmostoma branches with included Urceolus, and Peranema with the undescribed 'Jenningsia II', while other relationships are weakly supported and consequently the closest sister group to Euglenophyceae remains unresolved. Our results are inconsistent with recent inferences that Entosiphon is the evolutionarily pivotal sister either to other euglenids, or to Spirocuta. At least three transitions between posterior and anterior flagellar gliding occurred in euglenids, with the phylogenetic positions and directions of those transitions remaining ambiguous.

Phylogenetic and Morphological Evolution of Green Euglenophytes Based on 18S rRNA.

Green euglenophytes are a group of eukaryotes with ancient origin. In order to understand the evolution of the group, it is interesting to know which characteristics are more primitive. Here, a phylogenetic tree of green euglenophytes based on the 18S rRNA gene was constructed, and ancestral states were reconstructed based on eight morphological characters. This research clarifies the phylogenetic relationships of green euglenophytes and provides a basis for the study of the origin of these plants. The phylogenetic tree, which was constructed by Bayesian inference, revealed that: Eutreptia and Eutreptiella were sister groups and that Lepocinclis, Phacus, and Discoplastis were close relatives; Euglena, Cryptoglena, Monomorphina, and Colacium were closely related in addition to Trachelomonas and Strombomonas; and Euglena was not monophyletic. An ancestral reconstruction based on morphological characters revealed seven primitive character states: ductile surface, spirally striated, slightly narrowing or sharp elongated cauda, absence of a lorica, chloroplast lamellar, shield or large discoid, pyrenoid with sheath, and with many small paramylon grains. However, the ancestral state of the length of the flagellum could not be inferred. Euglena and Euglenaria, which both possessed all of the ancestral character states, might represent the most ancient lineages of green euglenophytes.